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ABSTRACT 

Assuming the weak energy condition, we study the nature of the non-central shell-focussing 
singularity which can form in the gravitational collapse of a spherical compact object in 
classical general relativity. We show that if the radial pressure is positive, the singularity 
is covered by a horizon. For negative radial pressures, the singularity will be covered if the 
ratio of pressure to the density is greater than —1/3 and naked if this ratio is < —1/3. 
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There are various known examples of formation of black-holes and naked singularities in 
spherical gravitational collapse (for recent reviews see [|I|]). However the general conditions 
on the initial data in spherical collapse which will ensure formation of a black-hole or of 
a naked singularity are not well-understood. In this Letter we demonstrate the role of a 
negative radial pressure in the formation of a naked singularity in the gravitational collapse 
of a spherical compact object. We use methods similar to those developed by Podurets @ 
and Lifshitz and Khalatnikov || to describe solutions of Einstein's equations in the vicinity 
of a singularity. 

In comoving coordinates (t, r, 8, <ft) the spherically symmetric metric is 

ds 2 = e a dt 2 - e w dr 2 - R 2 dQ 2 (1) 

where a, u and R are functions of t and r and R(t, r) is the area radius at time t of the fluid 
element with comoving label r. The energy-momentum tensor is = diag(—p,p r ,p0,pg) 
where p r and pg are respectively the radial and tangential pressures. We assume the weak 
energy condition: p > 0, p + p r > and p + pg > 0. The field equations are (see e.g. ||) 

m! = 4irpR 2 R' , (2) 

m = -Airp r R 2 R , (3) 

G H' 

2p e RRR' = -2m' + m'— + m— , (4) 

G H 

- 2R! + R'^ +R^ = 0, (5) 



m 

2 



]-R (l + e~ a R 2 - e^R' 2 ) . (6) 



Prime and dot stand for d/dr and d/dt respectively. We have defined G = e~ w R' 2 and 
H = e~ a R 2 . m(t,r) is an arbitrary function of integration, and from (fj) it follows, for a 
fixed time t, that 

m(t,r) =4rr J pR 2 dR + m(t,0). (7) 

m(t, r) is interpreted as the mass interior to r at a given time t. We assume that m(t, 0) = 
and that the total mass m(t,r Q ) is finite, where r is the boundary of the collapsing object. 
The weak energy condition ensures that m > 0. At the initial epoch t = t we assume 
p(ta,r) to be finite and strictly positive for r < ro and zero for r > tq. Eqn. ([?D along with 
m(to, 0) = then implies m(to, r) > for r > 0. 
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We assume that the gravitational collapse of the star leads to the formation of a shell- 
focussing curvature singularity. The singularity curve is given by R(t s (r),r) = 0, i.e. the 
shell with label r shrinks to zero area radius at comoving time t s (r). Further, we make 
the natural assumption that (dt s (r) / dR) ^ — different shells become singular at different 
comoving times t. Since we are considering collapse, we take R < 0, the equality sign 
holding, at best, initially and in the limit of approach to the singularity. The shell-focussing 
singularity at r = is called the central singularity, and that at r > non-central. In this 
Letter we are concerned with the non-central singularity. 

Assume first that the radial pressure p r remains non-negative throughout the star during 
evolution. It then follows from (J3|) and the initial conditions that m(t, r) is non-decreasing. 
For r > 0, m(t,r) is strictly positive. A shell labelled r, with r > 0, becomes trapped at a 
time tx{r) given by i?(i>r(r),r) = 2m(tT{r),r). Positivity of m implies that ty(r) < t s (r). 
The shell gets trapped before becoming singular and hence the singularity is covered. This 
argument does not hold for the r = singularity, since m(t,0) = 0. (In that case the 
ratio 2m/ R in the limit R — > 0, r — > will play a decisive role and in principle the central 
singularity could be naked.) 

Thus the non-central singularity will necessarily be covered if the radial pressure is pos- 
itive, irrespective of the sign of the tangential pressure. In particular, this includes the case 
of the radiative equation of state p r = Pe = p/3 which is expected to hold during the late 
stages of collapse. The weak energy condition, however, does allow for p r to be negative, 
in which case it follows from (|3]) that in the limit of approach to the singularity, m(t, r) for 
r > could go to zero and hence the non-central singularity could in principle be naked, 
unlike when p r is positive. This also means that any naked singularity arising in spherical 
collapse is necessarily massless, in the sense that the function m(t, r) will go to zero at the 
naked singularity. 

In order to illustrate how a negative radial pressure might alter the nature of the non- 
central singularity, we take the example of a perfect fluid with an equation of state p r = pe = 
p = kp, with — 1 < k < 0. We could assume that such a distribution might arise in the inner 
regions of a star, and is suitably matched to outer regions having positive pressures. Here we 
study only the implication of negative pressure, without considering how such a matching 
could be achieved. For a perfect fluid the field equations ([|) and (|5|) can be written as 

' 2? ' (8, 



P+ P 



2p 4R 

uo = . (9) 

P + P R v ; 

Eqns. (g), (Q) and (|) remain as before, with p r = p. (|8|) and (||) can be solved to yield 
e~ CT = p 2k / 1+k and = R 4 p 2 ^ 1+k , where we have used coordinate freedom to set constants 
of integration to unity. 
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Next we change from variables (t, r) to variables (R,r), and the three basic equations 
(0), fl3|) and (§) which have to be analysed further become 



^ = -^kpR 2 , (10) 

dt dm , , x dt , . 

__ = _ M i + ^, (id 

m ~2 R+ 2 Rp Jdt/Mf-2 \ Rp WJdRp- (12) 

Following Podurets || we make a crucial assumption, namely that the non-analytic part 
in p(R, r) depends only on R — this is physically motivated because there are no special values 
of r. The only exception is r = where this assumption is unlikely to hold; in particular it 
does not hold at r = for Tolman-Bondi dust collapse. Hence the following analysis is valid 
only for the non-central singularity, subject to the previous analyticity assumption. Thus, 
p = ip(R)pi(R,r) where p\ is analytic. By expanding px(R,r) in a power series near R = 
and keeping only the leading term we can approximately express p as 

p(R,r) = ip(R)a(r). (13) 

We assume that near the singularity curve we can write dt/dr fh dt s (r)/dr. Eqn. ( |T0"D then 
gives 

m(R, r) = -Airka(r)^(R) + m (r) (14) 
where d^/dR = ip(R)R 2 . Eqn. (O) can be written as 
dt dt s i)(R)R 2 

aR~ _7r( + 'li^ a[r ' {dm /dr) - Ank£(R)(da/dr)' { ' 

We first assume that as the singularity is approached, the mass function m(R, r) goes to 



a non-zero value. Then in ([12]) the second term should either be of the same order as the 
third term or dominate over it. Further, the second term should go to a non-zero value in 
this limit. It can be shown that if these two terms are of the same order then Eqns. (0), 
( |I4"D and ( |i~5| ) are inconsistent with each other. Hence the second term is non-zero in the 
limit, while the third term goes to zero. It can then be shown from flT2| ) and (|T4T ) that a 
necessary condition for the mass to go to a non-zero limit is that k > —1/3. If A; < —1/3 
the mass goes to zero in the limit. Note that for k > —1/3 the mass may or may not go to 
zero in the limit. If it goes to a non-zero limit, the singularity will be covered, as discussed 
above. 
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Consider next the case that the mass function does go to zero in the limit of approach 
to the singularity, and does so as m(R,r) ~ R n . From fllUD we get p ~ R n ~ 3 , and from (|ITD 
that (dt/dR) ~ R' 1 . Then (P) becomes 



m = Ifl + Ao ( r ) i2 3+2Kn-3)/(l +fc ) _ S ( r ) fl r+2(n-3)/( 1 +*) j (16) 

where A)( r ) and B (r) are positive quantities. In order for the density to blow up at the 
singularity, it is necessary that n < 3. This, along with the original assumption — 1 < k < 
implies that as R — > 0, the second term in (|16D is negligible compared to the first one, and 
hence in this limit (|16|) can be written as 



m=-R- B (r)R 7+2 ^ 3 ^ 1+k \ (17) 
2 

If in the limit R — > we have m = R/2, this is precisely the condition for formation 
of a trapped surface, and in this case we do not consider the singularity to be naked. It is 
at best marginally naked, in the sense that only one light ray escapes from the singularity. 
For the singularity to be naked, it is necessary and sufficient that in the limit i? — >• 0, 
(2m(R,r)/R) < 1. This requires that in (|17D the power of the second term should be equal 
to one, i.e. n = —3k. Further, for nakedness it is necessary that n > 1 (otherwise 2m/ R 
will be infinite in the limit R — > 0). This gives k < —1/3 as a necessary condition for 
nakedness. Hence for k > — 1/3 the singularity will be covered, irrespective of whether or 
not the mass goes to zero. 

For k < —1/3 the singularity is naked. This follows because if it were covered, then the 



power of the second term in (17) will be greater than one, giving n = 1 and hence k > —1/3. 
Hence we have the final result that for k < —1/3 the non-central singularity is naked, and 
for k > —1/3 it is covered. The result is subject to the crucial but plausible assumption 
about the analyticity properties of p(R, r) stated above. The role of negative pressure in 
naked singularity formation has been highlighted also in ||, although the starting point and 
some of the conclusions of that work are different from ours. 

It may appear that the condition 2m/ R < 1 need not be sufficient for naked singularity 
formation. For instance, it was suggested in |J that the absence of an apparent horizon 
up until the time of singularity formation indicates the singularity is naked. However, it 
was pointed out in [[7] and || that although one may not detect the apparent horizon at 
the time of singularity formation, in a particular space-time slicing, a trapped surface may 
nonetheless exist. This argument of and ||, however, is not applicable to our criterion 
above, because we are using the area radius R itself as one of the coordinates. Hence the 
result 2m/ R < 1 is independent of the slicing and indicates the absence of a trapped surface. 
Of course, the most direct way of showing the singularity to be naked is the existence of 
outgoing goeodesics. Nonetheless, our present criterion, though less direct, appears to be on 
the same footing. 
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One could ask whether a curvature singularity forms in collapse in spite of the negative 
pressure. For instance, it is well-known that in Robertson- Walker cosmology the evolution 
of the scale factor is given by 

S(t) = ~(p + 3p)S(t). (18) 

There is a past singularity S = if S < 0, i.e. p + 3p > 0, (i.e. k > —1/3 for p = kp). 
However a singularity may not form if p + 3p < 0, (i.e. k < —1/3). The weak energy 
condition (p > —p) by itself does not guarantee singularity formation in this case. However, 
consider the corresponding equation for R(t, r) in inhomogeneous collapse, which can be 
derived by differentiating (|) w.r.t. t and substituting from (|) and (§), (see also ||). 

R(t,r) = --Re' hp + J + - JL-^. (19) 



lq) can be recovered from (|19|) by setting i?(t, r) = rS(t), e a = 1 and r) = 1 — A;r 



( P~9]) differs from (18) in interesting ways, which suggest that a singularity could form even 
if p + 3p < 0. For densities decreasing outwards, the second term inside the bracket in (|19D 
is larger than p, allowing p to be less than — p/3 before R becomes positive. Similarly, the 
second and third terms in (|18"D are negative in the present case, again supporting singularity 
formation by keeping R negative. It is hence possible that a singularity can form for k < 
— 1/3, in which case it will be naked. On the other hand, if a singularity forms only for 



k > — 1/3 it will be covered. We note that there is actually a singularity theorem [I0| which 
assumes the weak energy condition. Lastly we mention that if our assumption of positivity 
of energy density is dropped, the mass function could in principle go to a negative value 
at the singularity. A negative mass region could arise for instance in the case of charged 
sphere solutions of the Einstein-Maxwell theory, leading to the Reissner-Nordstrom repulsion 
phenomenon 
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